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GENERALIZATIONS OF LAGRANGE'S EXPANSION COMBINED WITH LIGHTHILL'S 


TECHNIQUE FOR UNIFORMIZING SOLUTIONS OF 
PARTIAL DIFFERENTIAL EQUATIONS 
E. Dale Martin 
Ames Research Center 


SUMMARY 


Use of Lighthill's technique to uniformize approximate solutions of 
partial differential equations is simplified by incorporating it into a 
perturbation-expansion scheme based on a higher dimensional generalization of 
Lagrange's expansion. Lighthill's technique is made easier to apply by use of 
explicit formulas for the uniformized solution in terms of the previously 
determined nonuniform solution. Results of this study also indicate that 
Lighthill's technique can be more useful than realized previously. As one 
example, uniformly valid thin-airfoil solutions can be obtained, by the direct 
procedure, to any higher order of approximation. 


INTRODUCTION 


Since Lighthill introduced his technique for uniformizing approximate 
solutions of physical problems (ref. 1) , the method has been widely used, 
especially in fluid and gas dynamics, for both subsonic and supersonic flow 
problems. (See refs. 2-4 for discussion and references.) The basic princi- 
ples of the method have also been incorporated into extensions and other 
related approaches (refs. 3, 5, 6, 7). 

It has been widely believed for a number of years that the use of 
Lighthill's technique should be restricted to ordinary and hyperbolic differ- 
ential equations. No way has been found to apply Lighthill's technique in 
many problems governed by parabolic equations where the nonuniformity is due 
to an essential singularity, which can only be eliminated by "stretching" or 
magnifying the coordinate (cf. ref. 8). Difficulty in applying Lighthill's 
technique to solutions of elliptic equations has also caused the belief that 
its use was not generally valid there. Although Lighthill's technique had 
been used early to cope with the nonuniformity at the leading edge in thin- 
airfoil theory governed by an elliptic equation (ref. 9), it was subsequently 
presumed that Lighthill's treatment of thin airfoils was a special case, and 
that the method could not yield an improved uniform solution beyond the second 
order. Tsien (ref. 3) (among others) stated (incorrectly, as we shall see) 
that "a solution uniformly valid to all orders is not possible." This was 
cited as a "failure of the PLK method." Very recently, however, Hoogstraten 
(ref. 10) was able to use Lighthill's technique in conjunction with a special 
conformal-mapping technique to find uniformly valid thin-airfoil solutions for 



round-nosed and sharp-nosed airfoils. Essentially the same method has been 
used independently by Bollheimer and Weissinger (ref. 11) . These significant 
results obtained by Hoogstraten and by Bollheimer and Weissinger apply only to 
thin- airfoil theory governed by the Laplace equation. However, it will be 
shown that uniform higher order thin- airfoil solutions can be obtained by a 
divect application of Lighthill's technique in this elliptic problem. In the 
case considered, the solution is obtained quite easily in a variety of ways, 
and is possible to all orders. The procedure is not limited to Laplace's 
equation and so should be applicable to many other elliptic problems. It is 
applied in essentially the same manner to hyperbolic problems. (Conditions 
for applicability of Lighthill's technique are not investigated here. The 
procedures given are intended to apply only when Lighthill's technique is 
applicable. 2 ) 

It appears that the most simple and direct way of applying Lighthill's 
technique is to use a Lagrange expansion-perturbation scheme and to specify 
the terms according to Lighthill's principle that higher order solutions shall 
be no more singular than the first. The procedure has been given in refer- 
ence 13 for one independent variable. For extension of this approach to prob- 
lems where more than one independent variable must be strained simultaneously, 
a higher dimensional Lagrange expansion is needed. Several forms of the gen- 
eralization have been given in references 14, 15, and 16. Those necessarily 
involved derivations provide the general term and convergence criteria. Two 
independent derivations of simple procedures for obtaining explicitly the 
terms to any desired order of the N-dimensional Lagrange expansion in simplest 
form are given in chapter V of reference 17. The needed perturbation expan- 
sions based on Lagrange's expansion also have been given by Sack (ref. 18) in 
rather complicated form, in terms of multiple summations of a general term. A 
simple derivation of the needed terms to any desired order is given in 
reference 17 in a form most conveniently used for present purposes. 

Accordingly, the purposes of this paper are to outline a relatively 
simple and direct procedure for obtaining, to any order, the terms of a higher 
dimensional (vector) generalization of Lagrange's expansion; to give, in 
simplest form, a general perturbation-expansion scheme based on that generali- 
zation; to combine the expansion scheme with Lighthill's uniformization tech- 
nique for a considerable simplification of its use for any number of 
independent variables; and to illustrate the simplified application of 
Lighthill's technique in uniformizing the approximate solutions of partial 
differential equations. 

The special advantages of this approach to applying Lighthill's technique 
are: (a) By providing explicit formulas for the uniformized solution in terms 

of the original functions that were not uniformly valid, it eliminates much of 
the tedious procedure normally followed; and (b) it extends this simplifica- 
tion to higher dimensions, when more than one variable must be strained simul- 
taneously. For one independent variable, results of a previous paper (ref. 13) 
are extended to give a more general, and hence more flexible, expression of 

Comstock (ref. 12), in studying a problem posed by C. C. Lin, showed 
certain limitations to Lighthill's technique in ordinary differential equa- 
tions that presumably can apply as well to partial differential equations. 


2 



the explicit formulas for the uniformized solution. Because of this 
generalization, natural choices of constants in the uniformizing transforma- 
tion are often evident, and can lead more readily to the appropriate uniformi- 
zation. As a result of these factors it is believed that Lighthill's 
technique is made easier to apply, and more useful than previously, by the 
developments to follow. 

Helpful comments by Profs. M. J. Lighthill and R. A. Sack are gratefully 
acknowledged. 


VECTOR GENERALIZATION OF LAGRANGE'S EXPANSION 


Lagrange's Expansion in One Variable 

The standard form of Lagrange's expansion for one implicitly defined 
independent variable. 


Z = Z (?,e) - Z + ep(Z) 


(la) 


is 


fCZ) ■ fCO 

n=i 


d n - 2 



dfCQ I 
dc . 


(lb) 


Vector and Tensor Definitions and Notation 


For the treatment in higher dimensions, consider the N-dimensional space 
with orthogonal unit base vectors e k (k = 1, 2, . . . , N) : 


e i • e j = «ij = 1 
= 0 

Let , Z, and the function 
that 


for i 


(i » j = 1, 2, . . ., N) (2) 


for i t j j 

M(Z) be N-dimensional vectors in this space such 


Z = Z(J,e) = f + e/x( Z) 


( 3 ) 


where 

N 

? = E e k?k (4a) 

k=i 

• N 

Z = ^ 1 e k^k (^ i > ?? » * * ' * ^N 1 ^^ (4b) 

k=i 
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M(Z) = Xi e k^ J k^ 1> • • •» Z N) 

k=l 

For any arbitrary differentiable function F(?,e), define the vector 
operators and V z : 

N 


V E Fff,0 =]£e, 


Ftf ,E) 


V Z F(Z >E ) = s|-FCZ,e) 


where the partial derivatives are taken holding all other components of the 
argument fixed. 

For arbitrary N-dimensional vectors 


1 1 

A E ^2 e k A k ’ B E ]C 6 ^ 


we use an extension of the notion of an N-dimensional dyadic (second-order 
tensor) : 

N N 

AB = E E e i e 

i=l j=l 

to define the nth-order tensors 


A^ nJ e AAA . . . A , B 


(n) _ 


BBB 


n times 


n times 


We might call A*- n ^ and B^ "polyadics," since the special cases for n = 2, 
3, and 4 are known, respectively, as dyadics, triadics, and tetradics (cf. 
ref. 19) . The following defined scalar products then follow quite naturally 
from equation (2) : 


A • B e £ AiBi 
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N N 


AA; BB = A • (A • BB) = £ £ 

i=i j-i 


(9b) 


AAAiB BB = A • [A • (A • BB B) ] 

N N N 

= E E 53 A^AjAkB^Bj Bi 

i=l j=l k-i 


(9c) 


and, in general, define the nth scalar product: 


, , , . N N 

A (n) (?)B( n) = £ E 

i l = 1 i2 =1 


N 

E 

in =1 


A ii A i 2 


. . Ai B-! B 


in“ x n ul n-i 


ii x (9d) 


Particular examples of nth-order tensors to be used are: 

[m cr ) ] Cn:) s M(r)M(?) . . . mc?) 

' V 

n times 


N N 

E E 

i ! = 1 12 =1 


N 

E 

i n =i 


'ii 




p in cr) 


(10a) 


and 


,(n) 


- ^5 ' • * V S 


N N N 

■EE---E-*. 

i l = 1 i2 = 1 


an 


x n =1 


'!n 3?. ... 3?^ 

ii in 


(10b) 


N-Dimensional Lagrange Expansion 

The well-known N-dimensional Taylor expansion for a function f(Z) about 
the point ? is 


5 


9f(q,q,. 


f(Zl>Z2,. • • ,Zn) = 


N 


CZi - q) 




3q 


1=1 


N N 



i=l j=l 


3 2 f(q>q»* ' • >%) 

3^3^ 


(ID 


In the above-defined notation, with Z defined by (3), the complete 
expansion (11) is equivalent to simply 


f(Z) = fen ^cz)] (n) (?)vj n) f(n d2a) 

n=i 

Then, also, 

oo 

*t(Z) = M(?) + [M(Z)3 (n) (n)V^ n) P(?) (12b) 

n=i 


These two equations determine the higher dimensional Lagrange expansion to any 
desired order in e, as shown by the following example to order e 2 : 


f (Z) = f(f) + eju(Z) • V 5 f(f) + j e 2 [M(Z)M(Z)]:V c V c f(r) + 0(e 3 ) (13a) 


where 

M(Z) = mCT) + £M(Z) • V ?M (n + 0(e 2 ) 

= Mff) + eMtf) ‘ V^tf) + 0(e 2 ) 

Thus 

f(Z) = f(?) + eM( n * V ? f(f) + e 2 ([M(n • V ? M(f)] • Vgfff) 
+ y [M(f)M(f)]:V 5 V c f(f)} + 0(e 3 ) 


(13b) 


(13c) 
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To order e 3 , this standard form of Lagrange’s expansion extended to 
N-dimensions is 2 


f(Z) = f(?) + e{M • V^f} + 


e 2 

21 


(M ‘ V ? (M • v 5 f) + Cm • v c m) 


v ? f> 


+ 



CM • V C [M • V C (M • V 5 f) + (M • V^M) • V ? f] 


+ 2 (JJ. ■ V ? M) * V C (M * v ? f) + [M * V C (M * V ? M) 


+ (M • V ? M) • V ? M] • V 5 f} + OCe 4 ) (14) 

where the argument of each function on the right side of equation (14) is ? , 
and where Z is defined implicitly by equation (3) . To obtain the form (14) 
from the procedure illustrated in equations (13), several identities are 
needed. However, this is not essential, as the form (13) and higher orders of 
that form, with the definitions (9) and (10), are also readily usable. (An 
alternative, completely independent derivation of eq. (14) is given in 
ref. 17. The general terms are obtainable from the more complex derivations 
and results of ref. 16.) 


A PERTURBATION-EXPANSION SCHEME 


Any number of additional parameters may be included as arguments of 
f, M , and Z. Thus, for example, with only one additional parameter a, if 


z = Z(r, e ,a) s r + CM( Z,a) 


(15) 


then equations (12) may be replaced by 


f(Z,a) = f(?,a) [M(Z,a)J Cn) (?)vj: n) ftf ,a) (16a) 

n=i 


2 Note that in equations such as (13c) and (14) an operator does not act 

beyond the closing bracket of a pair inside which it is located. 
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where 


M(Z,a) = MCr,a) + ^5?- [M(Z,a)] n (")V ? (n) M(f ,a) (16b) 

n=i 


The more explicit form to order e 3 is the same as equation (14) except that 
the left side is f(Z,ot) and the functions M and f on the right side have 
arguments (?,a). 

Now suppose that f(Z,ct) and M(Z,a) can be expanded in power series in 
a: 


f(Z,a) = £ a k_1 f k (Z); M(Z,a) = £ a k "V(Z) 
k=l k=l 


(17) 


If these expressions are substituted into equations (16), or into the 
generalization of (14) , and if the results are specialized to the case where 
a = e, one obtains a useful perturbation-expansion scheme. To order e 2 , the 
results are: 


f (Z , e ) = fi(Z) + ef 2 (Z) + e 2 f 3 (Z) + 0(e 3 ) (18a) 


and 


f (Z , e ) = ficn + e{f 2 cn + m itn • v ? fi(n> 

+ e 2 jf 3 c?) + Mitn • v ? f 2 (f) + M 2 cn • v^cn 

+ (j)miCT) • V c [Mi(?) • V^l(?)] 

+ (j) [Ml O') • V ? Mi(n) • V ? fi(f)]j + 0(e 3 ) (18b) 

where 

z = ? + eMjtf) + £ 2 {M 2 (?) + Mi (?) • V^Mi (?) ) + 0(e 3 ) (18c) 


For the special case where ?, M, and Z have only one component each, £, y , 

Z, the results to all orders, derived from (1), are: 

oo 

f (Z,e) = e n ' 1 f n (Z) (19a) 

n= 1 
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and 


f(Z,t) 


£ e“ _1 «nCO 

n= l 




where 



(19b) 


(19c) 


As noted in the Introduction, equations (18) and (19) should be 
obtainable from the more complicated general forms in reference 18. These 
results have been derived here more simply and stated in these forms for con- 
venient use later. 


COMBINATION OF LAGRANGE EXPANSION METHOD WITH LIGHTHILL'S TECHNIQUE 


For a problem in which the terms of the approximate solution in the form 
of equation (18a) or (19a) have been found to be not uniformly valid, the 
basic principles of Lighthill's uniformization technique are: 

(1) To reformulate the expanded solution, transformed in terms of 
initially undetermined "strained coordinates," and 

(2) To specify the straining transformation in a manner that removes the 
nonuniformity, according to Lighthill's principle: 


Higher approximations shall be 
no more singular than the first. 


For an arbitrary straining transformation (arbitrary Mj(?) in eq. (18c)), 
any transformed expanded solution is given by equation (18bjto order e 2 , 
where equation (18a) is the expanded solution in the "unstrained" coordinates. 
Thus, part (1) of Lighthill's technique is done. It remains to specify the 
terms so that Lighthill's principle is satisfied. 

It is useful to consider first why the principle achieves uniformization. 
The expansion of the exact solution will have successively higher order 
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nonuniformities in the terms that are of higher order in z if the exact 
solution contains z in such a way that its expansion shifts a singularity. 
For example, suppose the exact solution contains g = (Z + e) -1 . This 
expression is singular at Z = -e. When e -»■ 0, the singularity moves to 
Z = 0. This is exhibited in the expansion about e = 0: 

g = g(Z + e) s (Z + e)" 1 = Z” 1 - eZ" 2 + e 2 Z" 3 + . . . (20) 

In a perturbation solution that would attempt to find g as a function of 
Z and e, such singular terms will appear. However, if one made a transforma- 
tion 

C = Z + e (21) 

the complete term representing g is simply 

g = gCO = r 1 (22) 

where 

Z = ? - e (23) 

with which g is uniformly valid at Z = 0. The transformation simply 
removed e from an expression (g) that is singular at a value of Z depend- 
ing on e. Thus, equation (20) contains e, but equation (22) does not, so 
that "solution" for g and Z in terms of ?, expanded in powers of e, is 
uniformly valid at Z = 0. With this observation, for the general problem, 
one then simply seeks a transformation that precludes the occurrence of the 
higher order nonuniformities (and therefore that eliminates the shifting of 
the singularity when the solution is expanded) by removing z from terms 
whose expansion would shift the singularity. This is accomplished by 
Lighthill's principle. It should be obvious then that Lighthill's technique 
should not be expected to remove the nonuniformity in a solution that is 
caused by an essential singularity , which is removable only by properly magni- 
fying the variables and not by a simple straining. Thus, in an extension of 
Lighthill's technique (see ref. 3), Kuo combined a "stretching" with a strain- 
ing of the coordinate in a boundary- layer problem. The nonuniformity caused 
by an (exponential) essential singularity occurs in problems where the highest 
order derivatives are multiplied by the small parameter, and so are lost in 
the perturbation solution (cf. refs. 2 and 8). 

In the classes of problems for which Lighthill's technique can be used, 
Lighthill's principle can be applied as follows: To first order as e ■+ 0, 

f(Z,e) = fi(?) in equation (18b). To satisfy the principle then, one can 
simply determine each n (?) so that the coefficient of each e n in equa- 
tion (18b) is equal to some constant a n times f i (§* ) . Thus, specify 

f(Z,e) = p(e)fi(n (24a) 
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where 


oo 

P(e) E 1 + X) a n e ” 


n=l 


(24b) 


The constants a n are chosen for convenience in each problem. They might all 
be taken to be zero, but it is just as easy to leave them unspecified a priori 
for possibly greater convenience. Thus, from equations (18b) and (24), one 
may evaluate the M n (?) to satisfy: 

f 2 ff) + Miff) • V 5 fi(f) = aififf) (25a) 

f 3 ff) + Miff) • v ? f 2 ff) + M 2 ff) • V^fiff) + (j) Miff) • VjMiff) • V 5 fiff)] 

+ (j) [Miff) * ^M iff)] • Vfiff) = a 2 f!ff) (25b) 


etc. 


For the reduction of equations (25) and (18) , the following notation is 
convenient. Let the component of Mff,e) in the direction of be 

yi(? 1 ,? 2 > • • ->?N> e )> which has the expansion 


ht (C ^ y • • *, ,e) l*il (s 1 » • • •> Cn) £ Fi2^1» • * b, ^n) * * * C^b) 

Thus, in y jj , subscript i indicates the component direction and subscript j 
indicates the order of approximation. One result is 

N 

Mj ff ) ’ ^5 - y ij [^1 > • • 9^ (h) 

i= 1 

for use in equations (25) and (18) . In the special cases where only one 
component of the independent variable (say, Zi) need be strained, denote: 


Hj 5 0 


— V*j(?l>« • •> £jj) 


for i f 1 
for i = 1 


If F* (? ) is then defined as 


F'ff) 


- 9F ff) 

Hi 


(28a) 


(28b) 
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we find, from equations (19b) and (24) , the special case of equation (25) when 
equation (28a) applies: 


Pi 

■Gr) 

( f 2 

- ajfi) 

( 28 c) 

P2 

II 

» —[ i — » 

['• 

- a 2 fi + Pif 2 + (j) Cuff{) 

( 28 d) 

^3 

'(if) 

f 4 

- a 3 fi + Pif 3 + P 2 f 2 + (j) ^l f z + 2y 1 y 2 f{)’ + ^ 

r) cu’fp"] 





( 28 e) 

^4 

= (*r) 

fs 

- a 4 f x + + y 2 fg + y 3 f 2 



♦(* 

)(P 

1^3 + 2 P 1 P 2 ^2 + 2l V 3 f ! + V 2 f ' 0 ' 



+ (k 

)<« 

1 f 2 + 3 PlP 2 £ P" + (^r) 

( 28 f) 


etc. In this special case, equation (18c) becomes (cf. eq. (19c)): 



The use of these equations for direct application of Lighthill's 
technique is illustrated below. With a nonuniform result in the form of equa- 
tion (18a) if more than one independent variable is to be strained, equa- 
tions (24), (25), (27), and (18c) are used; if only one coordinate is to be 
strained, equations (24) and (28) are convenient. 


APPLICATIONS IN IDEAL FLOW OVER A THIN BODY (ELLIPTIC EQUATION) 


The problem of flow over a thin airfoil of elliptical cross section is 
used for illustration here because it has been discussed extensively in the 
past in regard to the use of, and the previously presumed limitations of, 
Lighthill's technique. Also, the exact solution is available in a simple form 
for comparison of the results. 

The elliptical-airfoil surface is represented by 

y = y b (x) = ±e (1 - x 2 ) 1/2 , (-1 < x < 1) (29) 
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where x and y are inertial Cartesian coordinates, made dimensionless with 
respect to the length of the semimajor axis of the elliptical body. The flow 
is assumed to be steady, inviscid, and incompressible. The flow velocity is 
assumed to be tangent to the body surface and, far from the body, to approach 
a uniform stream velocity with magnitude U and direction parallel to y = 0. 

Consider the complex potential F(Z) and complex velocity W(Z) , which 
are analytic functions of the complex variable Z = x + iy: 


F(Z) = <Kx,y) + i<Kx,y) (30a) 

F*(Z) = W(Z) = u(x,y) - iv(x,y) | 

I (30b) 

= <t>x " 1< t‘y = 'J'y + i^ x ) 


where subscripts x and y indicate partial derivatives. The velocity poten- 
tial 4> and the stream function ^ then satisfy the elliptic equations 
representing conservation of mass and momentum: 


and 


^xx + 4* yy = 0 


$ ~ Ux + o(l) 


as x 2 + y 2 


(31a) 


x = iz. - 

u = * x 


dy b (x) 

dx 


on y = y b (x) 


v 

u 


'p 


XX + 'Pyy = 0 


ip ~ Uy + o (1) 
Px _ dyb(x) 
~ 4» v dx 


as x 2 + y 2 


(31b) 


or ij; = 0, on y = y fe (x) 


Consider the solution in three forms: F(Z), of which <p and $ are the 

real and imaginary parts; W(Z), of which u and -v are the real and imaginary 
parts ; and 


Qb(x) = [(u 2 + v 2 ) 1/2 l (32) 

L J y=y b U) 

the velocity magnitude on the body surface. Define the dimensionless 
perturbation-solution functions f(Z;e), w(Z;e), and q b (x;e) and the function 
g(x;e) by: 
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F(Z;e) = U [Z + ef(Z;e)] 
W(Z;e) = U[1 + ew (Z; e) ] 
QbCxie) E U[1 + eq b (x;e)] 

— U[1 + e + e 2 g(x;e)] 


(33a) 

(33b) 

(33c) 

(33d) 


A number of textbooks give the procedure for finding the approximate thin- 
airfoil solutions of equations (31). With y b (x) given by equation (29), one 
finds 


f (Z;e) = fi(Z) + ef 2 (Z) + s 2 f 3 (Z) + . . . ' 

w(z;e) s f’(Z;e) = f J (Z) + ef^(Z) + e 2 f 3 (Z) + . . . 

g(x;e) = g x (x) + eg 2 (x) + e 2 g 3 (x) + ... t 


(34) 


where 


fl(Z) = 

Z 

- (Z 2 

- D 1 ' 2 ; 

f 2 (Z) 

= fi(Z) j 

f 3 (Z) = 

fl 

(Z) - 

(!) t z2 - n " 1 ' 2 

; f4(z) 

= f 3 (Z)- ; 

f 5 (Z) = 

f 3 

(Z) + 

(!) » 2 - o’ 3 ' 2 

; etc. 


ref. 4) 

and 




g 1 (x) 

= 

-(f) 

x 2 ( 1 - X 2 )' 1 ; 

g 2 (x) = gj 

(x) ; 

g 3 (x) 

= 

(!) 

x 4 (l - x 2 ) -2 ; 

g 4 (x) = g 3 

(x) ; " 

g 5 (x) 

= 


) x 6 (1 - x 2 ) " 3 ; 

etc. 

i 


(35) 


(36) 


(cf. p. 52 of ref. 4) . 

For later comparison, note that the exact solutions (which can be found, 
e.g., from ref. 20, p. 429) are: 


f(Z;e) = (1 - e) ” 1 [Z - (Z 2 - 1 + e 2 ) 1/2 ] 


(37a) 
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(37b) 


w(Z;e) = (1 - e) ” 1 [1 - ZCZ 2 - 1 + e 2 )" 172 ] 
g(x;s) = — \ - 


We see that the approximate solutions (34) , with (35) and (36) , are 
singular at Z = ±1, the leading and trailing edges, and the higher approxima- 
tions are successively more singular at Z = ±1. We therefore consider use of 
Lighthill ! s technique by the direct procedure given above. To show the 
generality of the method, it is applied in the following paragraphs to this 
problem in a variety of ways. 


( r . x 

_\1 - x 2 + e 2 x 2 / 


(37c) 


First Solution 

If Lighthill f s technique in the form of equations (24) and (28), with the 
single independent variable Z, is applied to the complex potential F having 
the nonuniform solution given by (33a) with (34) and (35) , 

(e) (u" z ) e fCZ;e) = p(£)[e " a2 " 1)1/2] C38) 

If we now arbitrarily take all = 0 in equation (24b), so that p(e) = 1, 
equations (28c) to (28g) give the solution directly as 


f (Z;e) = C - U 2 - 1) 1/2 
Z = C + e(? 2 - 1) 1/2 


(39a) 


It is easily verified that this solution with only two terms in the expansion 
of Z is, in fact, equivalent to the exact solution (37a). Equations (39a) 
further reduce to 


^ = (1 + e)e 

Z = t, + £ (? 2 - 1) 1 7 2 


(39b) 


which has also been given by C. Jacob (see ref. 10). To find the expression 
for the corresponding velocity, one simply replaces f(Z;e) by w = f'(Z;e) 
and each f n (?) by f^(?) in equations (19). The result for w in (33b) is 
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w = 1 - S(c 2 - 1 )~ 1/2 + e[l - ?(? 2 - 1)" 1/2 - C? 2 - l)" 1 ] 


+ e 2 [1 - ?(c 2 - 1 r 1/2 - Cc 2 - 1) _1 ] + 0(e 3 ) (40) 


Note that there is no straining at Z = ±1, so the corresponding velocity 
would not be made uniformly valid by this straining transformation. Thus, 
even though this particular straining transformation (fixed by taking 
p(e) = 1) is successful in uniformizing the complex potential, it is not 
useful for the velocity. 


Second Solution 

If, instead of arbitrarily taking all a n = 0, one observes from 
equations (35) and (28c) through (28f) that it appears most natural to choose 
a n = 1 for all n (since it substantially reduces eqs . (28) for this problem), 
so that p(e) = (1 - e) -1 , one obtains directly from equations (24) and (28): 

f(Z;e) = (1 - e) _1 [C - U 2 - 1) 1/2 ] | 

z = 5 - ( x ) [5 " C ? 2 - l ) 1 ' 2 ]" 1 j (41) 

It is easily verified that this solution, with only two nonvanishing terms in 
the expansion of Z, is also equivalent to the exact solution (37a) . Further- 
more, there is a finite straining at Z = ±1. The corresponding velocity 
field, which then remains uniformly valid including the points Z = ±1, is 
given by (33b) , with 


w(Zje) = 1 - ?(? 2 - 1) _1/2 + e[l - Us 2 - 1)* 1/2 ] 

* s 2 [l - t(t 2 - l)- 122 - (i) (£ 2 - 1)-'] 

* e 3 [l - £t5 2 - l)" 1 ' 2 - (j) (C 2 - l)" 1 ] * OCE 1 *) (42) 

Third Solution 

If the velocity field is of prime interest, one might instead apply the 
uniformization directly to the complex velocity, W. The nonuniform solution 
is given by (33b) with (34), in terms of the derivatives of the functions 
listed in (35) . The uniformly valid solution can be represented by (24) and 
(28) with all f n (?) replaced by w n (?) = f^C?). Thus 

- l) = w(Z;e) = p(e)f{(?) (43) 
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where p(e) is given by (24b) and where the corresponding expansion of Z is 
given by (28g) . It appears most natural in this problem (in using (28c) 
through (28f) with all f n (0 replaced by w n (c)) to set a n = 1 for all n, 
as this substantially reduces the results. One then obtains directly 


w(Z;e) = (1 - e) ” 1 [1 - eU 2 - 1)' 1/2 ] (44a) 

Z - C - (~) eC - (|-) e 2 S - . . . (44b) 

It is easily shown, in fact, that the exact solution (37b) is equivalent to 
(44a) with Z = c(l - e 2 ) 1/2 , the expansion of which is (44b). 


Fourth Solution 

If one is interested in only the surface speed, (or the surface 
pressure coefficient, Cp = 1 - (Q^/ll) 2 ) as a function of x, one can use the 
above procedure directly to make uniformly valid. The nonuniform solu- 

tion is given by (33d) with (34) and (36) . The uniformized solution can then 
be found from equations (24) and (28) with Z and £ replaced, respectively, 
by x and 5, and with f, f n , and their derivatives replaced by g, g n , and 
their derivatives. Thus 

( 7 ^) [if " C 1 + £ )] E S( x ; e ) = P(<0g x (0 (45) 

where p(e) is given by (24b) and where the corresponding expansion for x as 
a function 5 and e is found from the form of (28g) . In using (28c) through 
(28f) one observes that it is convenient to let ai = 1 and a n = 0 for 
n >_ 2. One then obtains directly 

g(x;e) = (1 + e) ? 2 d “ 5 2 )" 1 ( 46a ) 

x = 5 + e 2 + e 4 (l - S 2 )" 1 + 0(e6) (46b) 


Not every possible choice of the straining transformation is useful in 
achieving the uniformization. The constants a n must be chosen to give 
useful results. 


Remarks on Elliptic Problems 

Because of the variety of ways the technique has been used (i.e., variety 
of straining transformations) in a straightforward manner in the perturbation 
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solution for potential flow over a thin elliptical airfoil, it appears that 
the technique should be applicable to a broad class of elliptic problems. 

Other simple examples of elliptic problems that have been worked out by the 
writer using the perturbation technique include potential flow over a thin 
parabolic cylinder, a thin round-nosed semi-infinite slab, a thin Rankine oval 
cylinder, and a slender axisymmetric paraboloid. Of course, it should not be 
implied that the procedure can apply directly to all elliptic perturbation 
problems. It cannot. The limitations of the method in elliptic problems are 
not definitely known, but are found to be apparently not so severe as believed 
previously. 


ILLUSTRATIVE APPLICATION TO A HYPERBOLIC EQUATION 


For illustration of the procedure when more than one independent variable 
must be strained to achieve a uniformly valid approximation, consider the 
problem: 


£ ** " £ yy + ( y£ x " x£ y) = °* < £ " °) 


(47a) 


y = 1, all x: f = 1 


(47b) 


3f _ £ — x 
3y “ e + x 


(47c) 


This problem is chosen for illustration because it is especially simple, 
and it combines several features that may be encountered in more complicated 
problems. Also, the approximate solution can be compared with the exact solu- 
tion. The partial differential equation (47a) is hyperbolic. Even though 
this particular problem could be simplified before solving it and the complete 
solution found, let us proceed to apply the perturbation technique directly, 
as one would need to do in more complicated problems. 

First, the perturbation solution is sought in the form: 


f(x,y,e) = f i (x,y) + ef 2 (x,y) + e 2 f 3 (x,y) + . . . (48) 

In the limit as e -*■ 0, one obtains the problem for fi, which can be written 
in the form 


3x 2 


CxyfO - £ 


(xyf i) = o 


(49a) 
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y = 1: fl = 1 


Hi 

3y 


-l 


(49b) 


The wave equation, (49a), with conditions (49b), has the solution 


fi(x,y) = y -1 


(50a) 


Similarly, the problem for f 2 is found to be: 



Oyf2) 


7-2 (xyf 2 ) 

9y 2 



y = 1: f2 = 0, 


af 2 

~w 


2 

X 


which has the solution 


f 2 ( x >y) = x-!(l - y~ 2 ) 


(50b) 


Further, one finds 


f3( x >y) = -x 2 y -1 Ci - y 2 ) 


(50c) 


The perturbation solution (48) with (50) is singular on both lines x = 0 
and y = 0, which are within the domain of interest. We therefore consider use 
of Lighthill f s technique according to the simplified approach outlined on 
pages 10 and 11. In particular, if the solution near the origin is of 
interest, both coordinates must be strained simultaneously. 

To apply this simplified procedure, let N = 2 in equations (4a) and 
(4b) and let 


z i = x, Z 2 = y 


(51) 


Then the expansion (48) corresponds to (18a) . The coordinate transformation 
(18c) implies (18b). To achieve the desired uniformization, equations (18) 
are also equivalent to (24a) with (24b) and with the functions [i . determined 
by (25). Thus 3 


f(x,y,e) 


(1 + aie + a 2 e 2 + 



(52) 
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where 


x — Zj — + ^ 


y = Z 2 = ^2 + ey 21 + £ l U 22 + 


(“12 “21 -5^)* ° ( “ 3) 

“11 T^- + “21 1^)* °t 


(53a) 


e 1 ) (53b) 


and where 


3f, 3f, 

f 2 * “11 ScT* “21 h 7- »l f l 


(54a) 


and 


/ 3£ 2 9£ 2 \ / 9f i 3f i \ 

f 3 + WT + U21 3^7/ + V Ul2 9 cT + y22 9CI/ 


1 / 3 3 \/ 3f i 9f i\ 

+ 2 ^31 3^ + V 2 1 a — 2 ) ^11 3^7 + *21 J^) 


+ 


1 

2 



3^11 

3?! 


+ U 21 


3 ’ J ll \ 3£i f 1. 

3? 2 / 3^2 ^ 



^21 + 
3?1 


P21 


3y 21 \ 9 ^1 
9? 2 / 3^2 


a 2 fi 


(54b) 


(the arguments of each fj and are c 1} £; 2 ). 

Equation (54a) gives an equation for y 21 in which there is no apparent 
advantage in having ai not zero; so for the simplest treatment, take 

ai = 0 (55) 


The result for use in (53b) 


If the determined functions 


is 

1^21 (?1>?2) 
(fi,f 2 , and 



1 


U 21 ) are then put into 


(54b) , 


(56) 
and if 


a 2 = 0 


(57) 
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for simplest treatment, the result for equation ( 53 b) to order e 2 is 


y = ? 2 + £ 


'4 - 1 ' 


+ 


^2 ~ ~ 1 
L \ , 


+ 0 (c 3 ) 


( 58 ) 


The function V*n(?i>?2) i- s as y et undetermined, but it can be conveniently 
specified from the following observation. We see that the expansion (58) for 
y 5 Z 2 is singular at Cj = 0 and that higher order terms in (58) would be 
more singular at = 0 if x : Zj were not strained appropriately. If we 

took p j 1 = 0, then Zi would be the same as Cj to order e, and the higher 
order singularity at = 0 in equation (58) would not be acceptable. We 
therefore invoke Lighthill's principle in a second sense in this problem, and 
specify to make higher order terms in equation (58) no more singular 

than (^2 - l)/?i at = 0. Thus, we may let 


?2 


~ ^11 


a = constant in equation (58) 


or 


V* 1 1 (? i » " “5] 


(59) 


and obtain the uniformized second-order solution: 


f(x,y,e) = — 

?2 


x = Ci + e(S 2 " + °C £;2 ) 


(60a) 

(60b) 


y = C 2 + e + 0(e 2 ) 


(60c) 


For the solution to be real very near the origin a must be less than zero; 
so we take 


a = ai = -1 


(61) 


Equations (60b) and (60c) are easily combined to obtain c 2 as a function of 
x, y, and e, so the uniform second-order solution for f is obtained directly 
as 3 

3 The appropriate algebraic sign on the radical in equation (62) or in 
(65) must be chosen, but this detail is considered incidental to the illustra- 
tion, as it is not of concern in equations (60) or (63). For both x >. 0 and 
y > 0, the + sign is used. 
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f(x,y,e) 


** |iLl_|y ± -L- [cx + ey)2 + 4e 2 (xy + e + e 2 )J | (62) 

If the above procedure is continued to order e 2 (leaving a as an 
arbitrary constant in equation (59)), the uniformized third-order solution is 
obtained: 


f(x,y,e) = 7 - (63a) 

x = + e (? 2 - ciCj) + e 2 a 2 ?j + 0 (e 3 ) (63b) 

+ 0(e 3 ) (63c) 

To this order, for the solution to be real very near the origin, a must be 
greater than zero; so we take 



a = a 2 = 1 (64) 

Equations (63b) and (63c) are easily combined to obtain so that the 

uniform third-order solution is expressed explicitly as 4 

f ** | * ~~ e v y ± jpr [( x + ey) 2 + 4e 4 (xy + e + e 4 )] J (65) 

For comparison, it is easily verified that the exact solution of the 
problem (47) is 


f (x,y,e) 


x + cy 
xy + e 


( 66 ) 


and that the approximate solutions (60) or (62), and (63) or (65), approach 
(66) uniformly as e -> 0 for all |x| and |y| < 00 . 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., August 12, 1970 


4 The appropriate algebraic sign on the radical in (62) or in (65) must be 
chosen, but this detail is considered incidental to the illustration, as it is 
not of concern in equations (60) or (63) . For both x >_ 0 and y >. 0 , the + 
sign is used. 
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